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The very strong electromagnetic fields present in ultraperipheral relativistic heavy 
ion collisions lead to important higher order effects of the electromagnetic interac- 
tion. These multiphoton exchange processes are studied using perturbation theory 
and the sudden or Glauber approximation. In many important cases, the multi- 
photon amplitudes factorize into independent single-photon amplitudes. These am- 
plitudes have a common impact parameter vector, which induces correlations be- 
tween the amplitudes. Impact-parameter dependent equivalent-photon spectra for 
simultaneous excitation are calculated, as well as, impact-parameter dependent 77- 
luminosities. Excitations, like the multiphonon giant dipole resonances, vector me- 
son production and multiple e'^e~ pair production can be treated analytically in a 
bosonic model, analogous to the emission of soft photons in QED. 



1 Introduction 

In ultraperipheral relativistic heavy ion collisions (UPC) very strong electro- 
magnetic fields are present for a very short time [1]. The parameter which 
characterizes the strength of the interaction is the Coulomb parameter t]: 



where Zi and Z2 are the charges of the nuclei and v ^ c is the relative veloc- 
ity of the ions. For heavy nuclei like Au-Au (RHIC) or Pb-Pb (LHC) r] ^ 1. 
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One can treat these systems semiclassically, with the nuclei following classi- 
cal Rutherford trajectories. At high energies they are well approximated by 
straight lines with impact parameter b. In addition to elastic scattering, which 
is due to the exchange of many photons, there are various kinds of inelas- 
tic processes, like nuclear excitations, photon-nucleus interactions leading to 
baryonic or mesonic excitations and photon-photon processes. One can calcu- 
late the amplitudes for these processes in the same way as the corresponding 
ones for real photons [2,3]. 

Probabilities for the excitation of the giant dipole resonance (GDR) and e+e"- 
pair production are very large, of the order of one for small impact parameter 
[2]. This means that these processes will also occur simultaneously with other 
interesting processes, where the probabilities are generally much smaller than 
one. The GDR is a strongly collective mode in nuclei. The theory of relativistic 
heavy ion excitation of one- and multiphonon GDR states is given in [4]. There 
is also mutual excitation, where each ion is excited to the (one phonon) GDR 
state. This process has been recently studied experimentally at RHIC [5]; it is 
important for the luminosity determination and for triggering on UPCs[6,7]. 
The fact that the double phonon giant dipole resonance is excited strongly in 
a two-photon excitation mechanism is of special interest for nuclear structure 
physics. The properties of these new collective modes could be explored with 
this method [8]. 

The probability of electron-positron pair production calculated in lowest order 
perturbation theory exceeds the unitarity limit of one [9,10,11]. This means 
that multiple e~^e~ pair production is expected to be appreciable for heavy 
ions at RHIC and LHC conditions. 

The probability P{b) of other inelastic processes is in general much smaller 
than one. An important process is diffractive vector meson production (espe- 
cially p°-production) . Probabilities for p"-production in close (grazing) colli- 
sions are about 0.01 or 0.03 for RHIC and the LHC, respectively. The proba- 
bility of the multiple (double) electromagnetic production of p° mesons is of 
the order of lO^^ to 10-^[12]. 

One example of simultaneous processes is giant dipole resonance excitation 
(followed by neutron decay) accompanied by p°-production. This process was 
recently studied experimentally at STAR[7]. 

This paper will present a general, unified treatment of multiphoton processes 
in heavy ion collisions, focusing on the correlations which occur because all 

of the processes share a common impact parameter, b. Section 2 reviews the 
semiclassical treatment and then introduces multiphoton interactions. Sec- 
tion 3 presents some examples of multiphoton processes and gives examples 
of the correlations, with RHIC and the LHC used for illustration. We give our 
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conclusions in Section 4. A preliminary account of part of the present work is 
given in Ref. [13]. 



2 Theory 



We work in the rest frame of one of the nuclei ("target system"). The other 
nucleus moves on a straight line with impact parameter b > Ri + R2, (where 
the nuclear radii are Ri and R2, respectively). The corresponding Lorentz- 
factor is 7 = is related to the Lorentz-factor ^coii in the 

collider frame by 7 = 27^^;; — 1. 

Since the strength of the electromagnetic interaction decreases with increas- 
ing impact parameter b, the excitation probability also decreases, in many 
important cases proportional to 1/6^. Higher order effects are proportional to 
the product of such probabilities, and are therefore especially important for 
collisions with impact parameter b close to — Ri + R2- 

The collision time is given by Tcou = b/'j, so the electromagnetic spectrum will 
include frequencies up to 

^max ~ ■ (2) 



2.1 Elastic Coulomb Collisions 



A simple and rather accurate result is well known for classical Coulomb scat- 
tering at small angles 9 [14] . The momentum transfer is 

2ZiZ2e^ 2r,h 



in the direction perpendicular to the motion. The momentum transfer is re- 
lated to the scattering angle 9 by Ap ^ p9. Prom this classical relation between 
the impact parameter and the scattering angle one can calculate the classical 
relativistic Rutherford scattering cross section as 



da 

In 



"2^1^26- 
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(4) 



where p = jMiv and Mi is the mass of nucleus Zi, see, e.g., p. 644 of [14]. 
This result can also be obtained in a field theoretical approach, which we only 
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Fig. 1. a) Classical picture of heavy ion scattering. This picture was also used by 
Fermi in his celebrated paper on the equivalent photon approximation [18]. He 
considered the nonrelativistic case u ^ c. The relativistic case is discussed in a very 
pedagogical manner in [14]. (b) For rj ^ I many photons are exchanged in elastic 
heavy ion scattering. A typical graph is shown. 

sketch here. Coulomb scattering is due to the exchange of photons. The sum 
over all the ladders diagrams, crossed and uncrossed, can be done in the high 
energy limit and for forward angles [15,16,17]. A typical diagram is shown in 
Fig. 1(b). 



One ingredient in this procedure is the linearization of the propagator denom- 
inator 



1 



1 



{p + q)"^ — w? 2p ■ q 



(5) 



In order to sum up all ladder and crossed ladder graphs one uses the identity 

1 1 



E 



^ 0(^(1) + a<T(2)) • • • («(t(1) + «(t(2) + \- 0'lO'2 " " " 



(6) 



where the sum is over all permutations a of the indices 1 to n. 



Summing over all orders the scattering amplitude is found similar to the 
eikonal approximation in (non relativistic) potential scattering: 



feiiq) = d%exp{iq- h) [exp(ix(&)) - 1] • 



(7) 



By integrating over the angle (f) 



fei{q) = ip J db bJo{qb) [1 - exp (ix(6))] 





(8) 



In this approach there is a well defined momentum transfer q = Ap (we 
set h — 1), and the impact parameter b is not an observable quantity. For 
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the elastic Glauber phase %(&) in the case of photon exchange, the analytic 
expression is well known, see e.g. [19]: 

X{b) = -2rjKo{Xb), (9) 



where A is a small photon mass needed to regularize the integral at large b. 

For 77 ^ 1 the integrand oscillates very rapidly and large contributions to the 
integral are only obtained when the phase of this integrand is stationary. One 
can introduce the saddle point ( "stationary phase" ) approximation [20] 



/ 27ri 



(10) 



where to is determined by $'(to) = 0. (One may have to sum over various 
stationary points to)- In our case we have a two-dimensional integration over 

— * 

6; for a derivation of the semiclassical case using only a one-dimensional saddle 
point approximation, see [21]. Assuming that the momentum transfer q points 
in the x-direction, the phase in the integrand of Eq. (7) is (for small A) 

$(6) = q^b^ + 2r]ln{Xb) (11) 



The conditions 

9$ ^ by b^ 



lead on the one hand to by = 0, that is the dominant contribution comes 

— * 

from the direction, where b is parallel to q. It also leads to an extremum for 
b^ — —bo with 

? = Ap=|^, 60 = ^. (13) 
bo q 



The second derivatives at the saddle point are 

-277 (92$ 277 , , 

= = = — (141 

dbl bl ' db.dby ' dbl bl ^ ^ 



Therefore an approximation of the elastic scattering amplitude up to a phase 
is 



pbl 2pri 

\fei[q) = ^ = — , 15 
2r) q'^ 
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giving the same cross section as the classical calculation. 

Eq. (13) gives the classical connection of the momentum transfer q to the 
impact parameter 60, see above Eq. (3). The dominant contribution to the 

— * 

(pseudo-) Gaussian integral comes from values of h around the classical impact 
parameter feo in (opposite) direction of with a width ("uncertainty of the 
6-value", see Eq. (14)) of ^. 

The classical momentum transfer Ap is built up from the perpendicular mo- 
menta of many exchanged photons. The perpendicular momentum of the in- 
dividual photon is restricted to ^ \-i ^'^ "^V photons must be exchanged 
to reach the classical Ap. This contrasts with electron or proton scattering, 
where rj = ^ 1 for v ^ c. So, in the following, we use the semiclassical 
approximation to describe inelastic processes in the case r] ^ 1. 

For b < Ri + R2 one has to consider strong interactions. In many cases, 
especially for heavy ions, the "black disk approximation" is quite good. In 
this approximation x is very large and imaginary for 6 < i?i -|- it!2. In the case 
of strong Coulomb clastic scattering (77 ^ 1), including the strong absorption 
of the ions for b < R1+R2, the diffraction pattern is given by Fresnel diffraction 
(rather than a Fraunhofer one for rj < 1). This is explained in detail in [22], see 
also Sec. 5.3.5 of [23]. These grazing angles are extremely small, about 5 x 10~^ 
(radians) for Au-Au collisions RHIC and 2 x 10~^ (radians) for Pb-Pb collisions 
at LHC (in the collider frame); they are well below the beam omittance, so the 
scattered ions will remain in the beam. Although the previous discussion is 
not entirely new, we think that it is useful to present the general ideas, which 
relate classical and quantal (or field-theoretical) descriptions of scattering in 
strong Coulomb fields, which may be a rather unfamiliar subject. 



2.2 Perturbation Theory for Inelastic Processes 

In the following we use the semiclassical approximation and treat the Coulomb 
field of the ion(s) as classical external fields. One ion moves on a straight line 
trajectory with a definite impact parameter, while the other ion is at rest. The 
semiclassical approximation can be justified in the same way as in the elastic 
case. The scattering amplitude in the inelastic case is given in the eikonal 
approximation by [21] 



00 




(16) 
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The amphtude o/j can be expanded in the form afi{b) = ajj(6) exp(i/x0) 
where fi denotes the angular momentum transfer along the beam axis. After 
integration over the azimuthal angle (p, 

oo 

IMq) = -W I dbY,MqbVexpMb)]a%{b). (17) 



Using the saddle point approximation, the connection to the semiclassical 
case can clearly be seen and a/j(6) corresponds to the semiclassical excitation 
amplitude. In the semiclassical approximation, the scattering amplitude is the 
product of the elastic Coulomb scattering amplitude and the semiclassical 
excitation amplitude afiib). We have 

ffM = feM)afi{-boq) = Uq) E(-l)'^<.(&o) (18) 



where q is the unit vector pointing in the g^direction (see Eqs. (11)-(14)). 
While we have not made numerical comparisons of the field theoretic solution 
to the semiclassical approximation in the relativistic case, such comparisons 
exist for the non-relativistic electromagnetic excitation. In this case the elec- 
tromagnetic cross sections can be calculated exactly. They depend essentially 
on the Coulomb parameter vj which approaches infinity in the semiclassical 
approximation. It is found (see, e.g.. Fig. (II. 9) of [24]) that the asymptotic 
(semiclassical) limit is reached already for rather low r] values. The semiclas- 
sical approximation works better than one may expect. In Au-Au or Pb-Pb 
collisions 77 ~ 50, so the semiclassical approximation should be excellent. 

The electric charge of the relativistic ion gives rise to an electromagnetic po- 
tential, the Lienard-Wiechert potential >l^(r, t) 



ylo(r, t) = 0(f, t) = ^^'^ 

^(b-xy + y^ + -f^(z-vty 

A{f,t) = ^<P{f,t). (19) 

The interaction with the target is described in terms of an electromagnetic cur- 
rent operator j^. This leads to a time-dependent electromagnetic interaction 
of the form 

Vit)^ I d'rA^{f,t)fif^. (20) 



The dependence of V{t) on the impact parameter is not shown explicitly. The 
target current describes a wealth of physics and is a quite complicated 
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object in general. This current contains the CTirrent of the nucleus and the 
interaction leads to the excitation of nuclear states, with the giant dipole 
resonance (GDR) as an important example. It also describes photonuclear 
interactions like nucleon excitations, meson production and photon-photon 
interactions, for example lepton pair production. The lepton (e'^'e") states are 
those in the external Coulomb field of the target, i.e., the Furry picture is used 
here. Now one can expand the (time dependent) target state in terms of 
all possible states which can be reached by the interaction V{t). We write 

^t) = ^ an{t) eM-iEnt)^n, (21) 

n 



where we have introduced the time-dependent amplitudes a„(t) for the state 
$„. The (time-independent) states consist of all states which can be con- 
nected to the target nucleus ground state by the interaction V (i.e., due to 
the interaction with the (virtual) photons). This includes, for example, excited 
nuclear states, the nucleus in its ground state and lepton pairs, or mesons, or 
nuclear excited states along with lepton pairs, or mesons, etc. 

One can set up coupled equations for the amplitudes a„(t) [25]. They are 
ian = Yl iM V{t) \m) exp(i(£;„ - Ejt)amit) (22) 

m 



The initial condition is am{t — oo) = Smo- The expression on the r.h.s. can 
be viewed as the matrix element of the operator 

V{t) = exp{iHot)V{t) exp{-iHot) (23) 



("interaction representation"). A formal solution of this equation can be writ- 
ten down using the time-ordering operator T: 

an(t ^ +oo) = (n| Texp j -z J dtV{t)\ \0) . (24) 



We use perturbation theory to solve this equation. Since the electromagnetic 
interaction V is in general weak, this is a good procedure. The first order 
amplitude is 

oo 

«n ^ ^ ~^ j dtVno{t) exp{iujnot) (25) 

— oo 
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where T4m(i) =< n\V{t)\m > and cunm — En — Em- This is the one-photon 
approximation. In second order 



1 ^ 



"n^^^Zl / ^^KmW exp j di'V^oit') e^p{iUmQt') , (26) 



— oo 



where the sum is over all intermediate states m. Only if the integrand is 
symmetric in t and t' one can extend the integral over t' to infinity and factorize 
the amplitude. This is not possible here since the operator V does not commute 
in general for different times t and t'. (In the case of the excitation of a 
harmonic oscillator, the commutator of these operators is a c-number and a 
full analytical solution can be found; this case will be treated below in Sec. 2.4.) 

In the case of the excitation of two independent modes, that is of two modes 
which do not mix with each other since their interaction is zero or negligible, 
there is such a factorization: let us assume that one excites a state n — (a, /3) 
where a denotes, e.g., a nuclear excited state and f3 a vector meson in a 
state with a given momentum and helicity produced coherently. Two types of 
intermediate states contribute to the sum over m, see Fig. 2: m = (0,/3) and 
m = {a, 0). By summing over them the integrand becomes symmetric in t and 
t'. Formally this can then be written in the following way: The first path goes 
from the state = (0, 0) to {a, 0) and then to {a, /?). The contribution of this 
path to the second order matrix element is 



^ oo t 

— j dtV(a,p)^a,o){t) ew {ii^(3ot) J dt'V^aMo,o)i't')^^Pii^aot') ■ (27) 

— oo — oo 

and for the second path 



«Sm(2) 



OO t 



^ J (iiV(a,;3)(o,/3)(t)exp(i(x;c,oi) J dt'V(o,/3)(o,o)(^') exp (io^/jo^') ■ (28) 



-oo — oo 



The independence of the two processes corresponds to the assumption that 
the presence of the state a (or /3) does not influence the interaction matrix 
element for the production of P (or a) (and that a single photon cannot 
produce a single step transition from (0, 0) to {a, j3), see below): 

V(a,/3)(0,/3) = V{a,omO) = Ko 
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GDR GDR GDR 



Fig. 2. Graphs contributing to the simultaneous excitation of states a and /3. In 
this figure we take a = GDR and = /9° as a typical example. 

and 

In many important practical cases this independence should be very well ful- 
filled. For mutual excitation the factorization is well understood, the elastic 
form factor of the excited ion is almost unchanged from the one of the ion in 
its ground state. Even the shght change in form factor takes time, as will be 
discussed below in connection with the validity of the sudden limit. 

Similarly for photon-photon processes like lepton-pair production together 
with nuclear excitation, the photon-photon process depends only on the charge 
of the nucleus, with little sensitivity to the details of the charge distribution. 
The nucleus is not influenced by the additional photon-photon process and vice 
versa. For hadronic flnal states, where there may be strong interactions with 
the nuclei, the photon-photon process in a double equivalent photon approach 
will predominantly occur outside the two nuclei. 

Adding the amplitudes of the two paths, Eqs. (27) and (28) one can see that 
the integrand is symmetric in t and t' and one obtains factorization. In an 
obvious way this is generalized to the excitation of N independent modes. For 
any function f{ti, t2,---, tn) which is symmetric in all variables ti 

t t 

j dti-- - j dtnfih, ■■■tn) ^ 
— oo — oo 

t tl *n-l 

n\ J dti J dt2--- J dtnf{ti,t2r--tn)- (29) 

— oo — oo — oo 

One can use this formula to factorize the amplitude into a product of flrst 
order amplitudes. The factor of ^ is compensated by the n\ different ways to 
order the excitation of the independent modes a, /3, 7, • • •. 
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(a) 



Fig. 3. Examples of graphs which spoil the factorization of the different processes, 
compare also with the figure above. As discussed in the text their contribution are 
expected to be small. 

Let us give here also an example of two modes, which are coupled (i.e. no 
longer independent): the excitation of vibrational and rotational modes of a 
nucleus [26] . These modes mix due to rotation- vibration coupling. This gives 
rise to nuclear eigenstates of a mixed character. In the case of p° production 
and GDR excitation, for example, such a coupling is expected to be so small 
that it can safely be neglected. 

The factorization is also not fulfilled in single photon exchange processes, 
which are depicted in Fig. 3. One such process is photon exchange with in- 
elastic vertices on both sides (Fig. 3(a)). This process was studied in [27,28] 
within an inelastic EPA approach and was found to be only of the order of 1% 
as compared to the elastic process. Another example is the GDR excitation to- 
gether with the meson production, see Fig. 3(b). This is a special contribution 
to the incoherent meson production process. Again this process is expected to 
be small: the vector meson is mainly produced through a diffractive (absorp- 
tive) process, which is very ineffective for a GDR excitation. The magnitude 
of such a process could be estimated by looking, e.g., at the GDR excitation 
in inelastic meson scattering at high energies. Unfortunately we are unaware 
of any study that has been done at these high energies. 

Asymmetric systems, such as d-A are also of interest, as shown by the recent 
d-Au run at RHIG, where UPGs were studied [29]. In these collisions, non- 
factorizable processes could come into play. The Coulomb parameter is 77 = 
0.6, less than one and the semiclassical approximation (used, i.e., in equations 
like Eqs. (3), (10) or (18)) does not apply. However, the total cross-sections (i.e. 
integrated over h and respectively) are equal in the semiclassical and Glauber 
(or plane wave) case [1,30]. In this run, the nuclear breakup channel was 
studied. Due to the dependence the photon(s) will come predominantly 
from the gold nucleus. The p° -|- breakup channel gets the main contribution 
from incoherent rho production, i.e. it is a one photon process of the form of 
Fig. 3(b). There is also a two-photon contribution from coherent p°-production 
accompanied by an electromagnetic deuteron breakup [31]. This mechanism 
is suppressed by a factor of Pbreakupib) relative to the incoherent one- photon 
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mechanism, which is of the order of 1% in the relevant impact parameter range. 
Due to the selection mechanism in the experiment the incoherent rho meson 
production mechanism is disfavored. Therefore overall the balance between 
the two reaction mechanisms should be studied in more detail. 



2.3 Multiphoton Processes in the Sudden or Glauber Approximation 



The solution of the coupled equations Eq. (22) is greatly facihtated if the 
sudden approximation can be applied, see, e.g., [30]. In this case the collision 
time tcoii = b/'-f is much smaller than the excitation time 1/uj. This condition 
is fulfilled in many interesting cases and we assume now that the sudden 
approximation can be used. This "frozen nucleus" -approximation is also used 
in Glauber theory. The relation between the semiclassical approach and the 
Glauber (or eikonal) approximation is explained in [30], for the non-relativistic, 
as well as, for the relativistic case. The excitation amplitude is 

an{t oo) =< n\ exp(zii:)|0 > (30) 
where 

+00 +0O 

R^- J dtVit) = - J dtA^{t)f. (31) 

— oo — oo 



The operator is a direct sum of operators in the space of nuclear states, the 
space of the nucleus-vector meson system, the nucleus-e"^e~ system, etc.. So, 

R = Y.Ri (32) 



where i denotes the different final states like e'^e pairs, p-mesons or a GDR 
excited ion. One has 

exp(ii?) = JJexp(ii?j). (33) 



as the different Ri commute with each other. One may expand the exponential 
in Eq. (30). Terms linear in R describe, e.g., the excitation of nuclear states, 
like the collective giant dipole resonance (GDR), vector meson production or 
e"'"e~pair production. Terms quadratic in R give, e.g., contributions to double- 
phonon GDR excitation, double vector meson production, double e'^e~ pair 
production. The quadratic terms also describe, e.g., vector meson production 
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and GDR excitation in a single collision. A contribution to the second order 
amplitude a^"^^ is, e.g., 



where |0 > denotes the ground state of the nucleus. In the general case the Ri 
commute and using Eq. (33) one gets 



which has the second order term as a special limit. 

For three independent processes, say GDR excitation, vector meson- and e+e" 
pair production, the six different time orderings compensate for the 1/3! factor 
in the expansion of Eq. (30). In this formalism, processes are independent and 
the elementary amplitudes factorize, as one would have intuitively expected. 
This property is used also in the experimental analysis of vector meson pro- 
duction with simultaneous GDR excitation, where the neutrons from the GDR 
decay serve as a trigger on UPC [7]. The ion motion is not disturbed by the ex- 
citation process. The reason is that the kinetic energy of the ion is much larger 
than the excitation energy. Due to coherence, the quantity x — ^ — ^^J^^ is 
very small. One has a; < 4 x 10~^, 3 x 10~^ and 1.4 x 10~^ for oxygen, tin and 
lead ions respectively. 

The GDR excitation immediately affects the nucleon momenta, but the nucle- 
ons will take a finite time to adjust their positions. With an excitation energy 
of about 20 MeV, the oscillation time of the GDR can be estimated classically 
to be about 10 fm/c. Due to the time dilatation, this corresponds to a time in 
the collider frame of about 1000 fm/c at RHIC and even longer at the LHC. 
This time is long compared to the time required to, for example, produce a 



From the general discussion above of higher order theory we have seen that 
this is not a necessary assumption to achieve factorization. 

2.4 Excitation of a Harmonic Oscillator (Bosonic Modes) 

An especially simple and important case is the excitation of a harmonic os- 
cillator. It also has the virtue that it can be studied analytically. As some 
excitation processes in heavy ion collisions can be regarded as quasibosonic 
processes, such a model can serve as a good approximation. Some examples 
are discussed below. Although well known (see, e.g., [32,33]), we briefly present 



a(2) = _ < GDR\R\Q >< p°\R\0 > 



(34) 



a =< GDR\ exp{iRGDR)\0 >< P°| exp{iRpo)\0 >, 



(35) 
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the main results of the excitation of a harmonic oscillator by an external force, 
this serves also to establish the notation. 



In terms of the creation and destruction operators and a the Hamiltonian 
of the system is 

H ^uj(^a''a+^^ (36) 



where cu is the energy of the oscillator. The boson commutation rules are 
[a,a^] — 1 and [a, a] — [a^,a^] — 0. Only one mode is shown explicitly, in 
general one has to sum (integrate) over all the possible modes. 

One assumes that the operator V (see Eq. (20) is linear in the creation and 
destruction operators: 

V{t) = f{t)a + r{t)al (37) 



In this case one can calculate V explicitly using the boson commutation rules 
given above and the expansion 

V{t) = V{t) + zt[Ho, V{t)] + ^[^0, [Ho, V{t)]] + ... (38) 



One finds 

V{t) = f{t)e-'^'a + f*{t)e"''a^ (39) 



Now one can convince oneself that the commutator of V at different times t and 
t' is a pure c-number. In this case one can disregard the time ordering operator 
in Eq. (24) and obtain an exact analytical answer, up to an unimportant overall 
phase factor (see, e.g., [32]). 

The excitation operator is the integral of the interaction over time 

+ 00 

J dtV ^ -{u a + u* a^) (40) 



where u is the c-number 

oo 

u= j dtf{t)exp{-iujt). (41) 
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This leads to the excitation of a so-called coherent state, see [34,35]. For the 
excitation of multiphonon states this is explicitly shown in [36]. One has 



an =< n|e-'("*"^+"«)|0 >= ^^^e'^"''*, (42) 



where the operator identity e^"*"^ = e^e^e"^!^'-^] was used, which is vahd 
for two operators A (— i?7,*a^) and B [—iua) for which the commutator is a 
c-number. See also the discussion of the forced linear harmonic oscillator in 
[33]. 

Let us look at a few cases where this model can be applied: Electromagnetic 
excitation of nuclear states, especially the collective giant multipole resonances 
was discussed recently by Bertulani [37]. The possibility to excite multiphonon 
GDR states is studied in [38]. Multiphonon GDR also play a role in the electro- 
magnetic excitation of the ions in relativistic heavy ion collisions, for a recent 
reference see [39,40]. The parameter which describes the probability PaDnib) 
of GDR excitation is (see [4]) 

m = I (43) 



with a simple parameterization for S for the relevant high beam energies [1] 

S = ^"'^'^"^^^ = 5.45 X 10-' ZlmZ.Af/'W. (44) 
A2mNUJ 



where itin denotes the nucleon mass, and the neutron-, proton-, and mass- 
number of the excited nucleus are N2, Z2, and A2 respectively. The excitation 
probability is inversely proportional to the energy cu (pa 80MeV/l~^/^) of the 
GDR state. As expected, the soft modes are more easily excited. Eq. (43) is 
based on the assumption that the classical dipole sum rule (Thomas-Reiche- 
Kuhn sum rule) is exhausted to 100%. For the excitation of an A^-phonon 
state, a Poisson distribution is obtained. For heavy systems P{b) ~ | for close 
collisions {b ^ Ri + R2). 

Double p° production was studied in [12]. In addition to the label m for the 
magnetic substates of the GDR, there is a continuous label (the momenta). 
The probability to produce a p°-meson in a close collision is of the order 1-3% 
for the heavy systems. It will be interesting to study these events experimen- 
tally to see how the p — p interaction affects the simple harmonic oscillator 
description, especially for the production of p°'s with close enough momenta. 

Multiple e"'"e~-pairs can be produced in relativistic heavy ion collisions [41]. 
This work used the sudden (or Glauber) approximation and a quasiboson ap- 
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proximation for e^e~ pairs. Thus the Hamiltonian has the form of Eq. (36), 
where a sum over the quantum numbers of the lepton pair has to be in- 
cluded. Using a QED calculation (including Coulomb corrections in the Bethe- 
Maximon approach) for one pair production as an input, a Poisson distribu- 
tion is obtained for multiple pair production. This is quite natural, since this 
problem reduces to the excitation of a harmonic oscillator (the modes are la- 
belled by the spins and momenta of the e+e" pairs). More detailed calculations 
have verified that the Poisson distribution holds quite well for multiple pair 
production, although shght deviations are expected [10,2]. A characteristic di- 
mensionless parameter for this problem is S = [{ZiZ'2.a^Y]/{mehY , where mg 
is the electron mass and h > l/mg. Here, log 7 should be large. 

With heavy systems like Au-Au or Pb-Pb, when h ^ 1/me, 2^1. The impact 
parameter dependence of e'^e~ pair production has been studied numerically 
in lowest order QED [9]. Only recently an approximate analytical formula for 
the total pair production probability in lowest order P^^^ was found. When 
1/me < 6 < 7/me [42] 

28 

P^^) = 7rT^(2 In 7^ - 3 ln(me6)) ln(me6) (45) 



The A^-pair production probability decreases strongly with increasing impact 
parameter h (approximately like ~ 6"^-^, for h > l/nie). Therefore the prob- 
abilities P^^\h) should be known accurately for an impact parameter range 
of < 6 < several l/mg. Since the multiple pair production happens in a 
single collision, with a given impact parameter, there are correlations in the 
momenta of the outgoing pairs. This will be briefly discussed in Sec. 3 below. 

e~^e~ pair production is of great practical interest. It could be useful as a 
trigger for UPC collisions at the LHC [43], but also constitutes a background 
for the most central parts of the trigger system of ALICE. [44]. Electron or 
muon pair production might also be a way to measure the 77-luminosity [45] . 

Muon pair production is also of interest: in close collisions the pair production 
probability is appreciable, although significantly smaller than one. This can 
be seen from Eq. (45). The muon Compton wavelength 1/m^ is 1.86 fm, much 
smaller than the nuclear radius. For b ^ 2Ra form factor effects also need to 
be considered. Total cross sections are smaller by about a factor of (mg/m^)^ 
^ 2x 10~^ (which is still large). The results of a lowest order external field QED 
calculation are shown in Fig. 4. Calculations with and without a monopole 
form factor are shown together with the approximate expression of Eq. (54) 
below. Details of the calculation are given in [10] . 
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Fig. 4. The probability of muon pair production is shown as a function of impact 
parameter for Pb-Pb colUsions at the LHC (7 = 3400). Results for a monopole 
form factor F{Q) = 1/(1 + Q^/A^) with A = 83 MeV and without a form factor 
{F{Q) = 1) are given. 

3 Some Examples of Multi-Photon Processes 

The factorization of the amplitudes has some significant, experimentally ac- 
cessible implications. Although the amplitudes arc independent, the processes 
all share a common impact parameter b. This leads to some significant cor- 
relations, due to both the magnitude and direction of b. In this section, we 
consider the effect of multiple interactions on the photon spectrum and impact 
parameter distribution, and on the angular distributions of the final states. 



3.1 Impact parameter dependence of the equivalent photon spectrum 
The photon spectrum at impact parameter b > Ra is 



where x = ujbj^. As long as x ^ 1, the photon spectrum scales as dN/dw ~ 
1/b^. The mean impact parameter b is 




(46) 



6 = 



fd% bP{b) 
Jd^Pib) 



(47) 
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where the probabihty P{b) of a specific reaction occurring at impact parameter 
b is 



Pib)^N{u;,b)a^A. (48) 

For a reaction involving a single photon with energy uj, we can integrate over 
the region 6^j„ = '2Ra < \b\ < bmax — ll^- Many processes can occur over a 
range of photon energies; in this case, an integral over to is needed. However, 
the details of the interaction do not matter, and 



For single photon processes the average impact parameter b is large when bmax 
is large and insensitive to the exact value of bmin- 

However, for two or more photon interactions, such as mutual Coulomb ex- 
citation of two nuclei, P{b) a; {(PN^/d%Y ~ 1/^^", with n the number of 
photon interactions and the latter approximation requiring that the x for the 
two photons both be less than one. Then, as long as bmax ^ ^mm, 

2n — 6 



bmax drops out, and the median impact parameter only depends on bmin, inde- 
pendent of the collision energy and other details of the interaction. Of course, 
bmax depends on the photon energy and we assume that for all photon pro- 
cesses bmax 3> bmin should be fulfilled in the relevant uj range. For heavy nuclei 
hke Au 62 ~ '^bmin — 4-Ra = 26 fm. 

For three-photon interactions, such as vector meson production accompanied 
by mutual Coulomb excitation of both nuclei, n = 3 and 63 ~ 46^i„/3 = 
For heavy nuclei, 63 ^ 18 fm. Of course, as 6„ drops with increasing n, the 
cross sections for UPCs accompanied by hadronic interactions become more 
sensitive to the tails of the nuclear density profile — including these tails 
increases bmin above the canonical 2Ra- 

The b are comparable to the median impact parameters calculated in [46]; 
these detailed calculations found that the median impact parameter is almost 
independent of the specifics of the reaction. This impact parameter distribu- 
tion is important for studying interference between vector meson production 
on the two nuclei; the smaller 6, the easier it is to observe the interference [47]. 

The reduction in b affects the photon spectrum, since the maximum photon 
energy scales as 1/6. One way to quantify this is to consider a two-photon 
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reaction: excitation of one nucleus to a GDR, while the second photon produces 
a p meson. The GDR tagging has a strong effect on the p'^-production: it selects 
small impact parameters, where the equivalent photon spectra are harder. One 
can calculate the equivalent photon spectrum in the case of GDR tagging. For 
close collisions, which matter most, the probability of GDR excitation is given 
in Eqs. (43) and (44). The equivalent photon spectrum is then 

Tl '^V J X \ 1 J 



This expression can be evaluated numerically. Due to the high values of the 
Lorentz boost 7 at the colliders, Kq{xY /'^'^ is very small and can safely be 
neglected. A simple approximation [48] is found by setting x'^Ki^xY — 1 
for X < 1 and zero otherwise; the integral over x is then —4 \ where 



Xmvn = < 1- Then 



nWUu.)^S^-fi^f[^-l\. (52) 



The photon spectrum scales very roughly as nioj) ~ a;°, in contrast to the 
n(a;) ~ a;~^ for untagged photons. 

Using Eq. (11.3.30) of Ref. [49] one finds a better analytic expression 

S Z'^ot 

riGDRiuj) ~ -^2 —^mini^li^minf " Ko{Xminf)- (53) 



Figure 5 shows the EPA spectrum n{Lo) (without tag) and the spectrum 
nGDRi^^) in the approximations of Eqs. (52) and (53) respectively. 

Another example of multiphoton processes is GDR tagging of 77-processes: if 
the mass of the system produced in the 77 interaction is low, there is a con- 
siderable influence of the tagging, if the mass is high, the influence is less. The 
6-dependence of (untagged) photon-photon-processes can be found from the 
folding of 6-dependent equivalent photon spectra, see e.g. Eq. (49) of Ref. [2]. 
A simple approximate formula was developed recently for the special case of 
e~^e~ pair production. The result is given in Eqs. (0.23), (0.16) and (0.17) of 
[42]. 

We consider here the more general case where the invariant mass of the 77- 
system is given by W (rather than 2me, as in Ref. [42]). The lower cutoff in the 
impact parameter bmin is 2i?^, (rather than the Compton wavelength ^ of the 
electron). We introduce the parameters u — ^ and x — . We reconsider 
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Fig. 5. The equivalent photon spectrum n{uj) (without tagging) is compared to the 
one for the excitation of a GDR resonance in one of the ions for Au-Au {Z = 79, 
A = 197) coUisions at RHIC {'^coll = 108). The full expression is compared to the 
two approximations given in Eqs. (52) and (53) 

Eq. (0.24) of [42], with different integration limits: bi from Rto b and Ui from 
W^R/A-f to 7/i?. Generalizing Eqs. (0.16) and (0.17) of [42], the ^^dependent 
77-luminosity is 

d^L^j _ 27r f Z'^a 
d?bdW ~ WR? 




The probability PAA{b) of a 77-process in the AA collision is given in terms 



of this 77-luminosity by PAA{b) 
the elementary 77 cross-section and 



"^^"J W d?bdW' 



(7-y-y. In this expression a^^ is 



X 



ln(a;^/'u) 



iOT X < u < 



(55) 



(for "large "impact parameters) . For "small" impact parameters 

X ^2l-a{u){2hi{x/u) + l/2hi{u)) ior 2 < u < x. (56) 

This result is obtained by a generahzation of Eq. (0.24) of [42]. 

The 6-dependence is strongly influenced by the value of W: for WR <^ 7 (e.g., 
for low mass e+e" production) there is an approximate 1/6^- dependence over 
a considerable range of 6-values, for WR ~ 7 the value of x is quite small 
and the dependence is much steeper. In this case low 6-values are emphasized 
anyhow and the effect of tagging is less dramatic than for the case of p"- 
production with GDR tagging, as discussed above. This was found by Klein 
in [50] for photon-photon processes with nuclear breakup. In accordance with 
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Fig. 6. The photon-photon luminosity d^L/d?bdW is shown for different values of 
W as a function of the impact parameter b. The results of a full calculation for lead 
beams at the LHC are compared with the approximate result of Eq. (54). 

the present considerations, the effect of tagging is more important for low W, 
see his plot as a function of the invariant mass of a produced muon pair iW^^. 
A comparison of the luminosity of Eq. (54) against a full calculation is shown 
in Fig. 6 for different invariant masses W. In view of the crude approximations 
made to obtain Eq. (54) the overall agreement is quite good for small invariant 
masses, whereas for the highest invariant masses shown in the plot some part 
of the cross section come from the areas beyond the range u < x'^. This figure 
also shows the change in the shape of the spectrum for low and high invariant 
masses W. 

3.2 Angular correlations 

There can also be angular correlations since the photon polarization follows 
the electric field vector of the nuclear fields which in turn follows b [50] . For 
example, the decay po ^ n~^7r~ is sensitive to the photon polarization, and 
hence to the direction of b. Mutual GDR excitation is another example; the 
neutron transverse momentum tends to follow the electric field. The amplitude 
ai,2(&) for mutual excitation is 



where the amplitude for single GDR excitation (i = 1, 2) is given, e.g., by 
Eq. (3.1.22) of [1]. 



In intermediate energy heavy ion scattering (in contrast to the ultrarelativistic 
case) it is experimentally possible to measure directly the momentum trans- 

— * 

ferred to the projectile. This determines the impact parameter vector b. The 



ai,2(b) = ai{b) 02(6) 
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decay products show a dependence on the azimuthal angle 0. This dependence 
has been observed, e.g., in the breakup reaction ^^Be + Pb -^^^ Be + n + Pb 
in [51], see their Fig. 3. The angular distribution follows a dipole distribution, 
see Eq. (3.1.22) in [1]. For large values of 7 the m — ±1 components dominate 
and 

Qi ~ sin 9 cos (f), (57) 



where 6 and are the polar and azimuthal angles of the emitted neutrons in 
the system of the nucleus. 

— * 

For mutual GDR excitation, the direction of b is the same for both excitations, 
producing a correlation between the emission angles of ^1,^2,^1, and 02 (in 
the system of the respective emitting nucleus). The correlation in the relative 
azimuthal angle A0 = 0i — 02 of the emitted neutrons is 

C(A0) = 1 + ^cos(2A0). (58) 



This correlation is shown in Fig. 7. The neutrons from the GDR decays can 
be detected with zero degree calorimeters which are located downstream of 
the collision region at relativistic ion coUiders. At RHIC, the ZDC's are lo- 
cated ±18 meters downstream; they are 10 cm wide by 18.75 cm high [52]. 
The dimensions are limited by the available space. Although the sensitivity 
is reduced near the edges due to transverse shower leakage, coverage is good 
for angles 9 < 2 mrad [5]. The excitation energy for a heavy nucleus GDR is 
approximately 10 MeV, corresponding to a neutron momentum of about 140 
MeV/c. 

For the RHIC neutron longitudinal momentum of 100 GeV, pt = 140 MeV/c 
corresponds to a deflection of 2.5 cm, well within the ZDC acceptance. With 
position sensitive shower-maximum detectors to the ZDC's, neutron position 
resolution of order 1 mm should be achievable [53], giving adequate angular 
resolution to study neutron-neutron correlations between the ZDCs. At the 
LHC, the ZDCs are located ±116 meters from the collision regions, giving a 
maximum neutron deflection of 5 mm, so determining the neutron angular 
direction may be more difficult. 

Similar angular correlations will occur for p production and decay to tt+tt" 
accompanied by single or mutual GDR excitation. In its rest frame the p decay 
angular distribution is given by \Yim{9, $)|^, where M denotes the magnetic 
quantum number. The population of the different M-substates depends on the 
dynamics of the production process. 

In the heavy ion case the population is relatively simple, since the photons are 
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Fig. 7. The angular correlation of two neutrons emitted after mutual GDR excita- 
tion is shown. denotes the azimuthal angle between the two neutrons 

almost real, with both helicities equally populated. It is the simplest assump- 
tion that in p photoproduction the p has the same polarization as the initial 
state photon - this is s-channel helicity conservation (SCHC). This has been 
experimentally verified for photoproduction [54,55] and for clcctroproduction 
[56]. In our case the photons are linearly polarized, transverse to the beam. 
Assuming SCHC, the final p° state is also linearly polarized, in the same di- 
rection as the equivalent photon. It is a superposition of M = ±1 states. The 
angular distribution in the rest system is 



quite analogous to Eq. (57) for the GDR decay. For exclusive p production, the 
direction of the impact parameter is unknown, and one has to integrate over all 
impact parameters. However, if the coherent p-production is accompanied by a 
GDR excitation, the azimuthal angles of the neutron and the Tr^vr^-pair would 
be correlated, as with mutual GDR excitation. The GDR decay provides a way 
to measure the photon polarization. The same correlation would be observed 
between two p° that are produced together [12]. 

In e'^e" pair production one also expects a dependence on the azimuthal 
angle which is defined by the impact parameter vector b. For polarized 
(real) photons this can be seen from the analysis of, e.g., [57]. For linearly 
polarized photons the electron is most likely to be emitted in the plane of 
the polarization [57]. It will be of interest to do such theoretical calculations 
in more detail and work on this is currently in progress. E.g. in two-pair 
production one might expect to see such correlations, which leave their trace 
on the impact parameter of the collision. The formula for one-pair production 




sin 6 cos $ 



i sin 6 sin $ 



(59) 
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in lowest order semiclassical approximation is given e.g. in [25,10,58]. This 
formula contains the correlations of the impact parameter and the momenta 
of the outgoing e+e~ pair. 

We conclude this section with a general remark on multiphoton processes. The 

observation of two electromagnetic processes with photon energies uoi and UO2 
gives information about the range of impact parameters of the collision. One 
has 

Pi2{h) = Pi{b)P2{b) = Niui, h)N{u2, h)a^{u^)a^{u2) (60) 

We use a simple approximation for the 6-dependent equivalent photon spec- 
trum N{u!,b): 

7?ni 

= for2i?<6<^ (61) 

Thus such an interaction takes place in the b-range of 2R < b < — where 
^max is the larger of uj\ and uj2- 



4 Conclusions 

Due to the strong fields of highly charged ions multiphoton processes are 
abundant in UPCs. The important modes were identified. Quite generally, the 
excitation of the different modes occurs independently and the corresponding 
amphtudes can be factorized. They do not influence each other. This is es- 
sentially due to the very high energy of the heavy ions, which is much higher 
than the energy uj of the equivalent photon. Only a small fraction of the total 
kinetic energy of the nuclei is lost in the UPC. Due to the coherence condition 
one has x^ax — ^ — rmj^a ~ ^ nucleon mass, and Ac (A) 

is the Compton wavelength of the nucleus. 

The main result of this paper may seem intuitively obvious. However, the ex- 
perimental implications are significant. The factorization assumption was used 
in detailed theoretical calculations of p-production and GDR excitation in [59] 
(following Table 2 it is said there that "the validity of the assumption of fac- 
torization is hard to prove rigorously ..."), and [46] and in the corresponding 
experimental analysis [7]. It is important to show also how this factorization 
arises in the theoretical formulation of the process. Now there is a formal basis 
to discuss these questions theoretically. 

Ultraperipheral heavy ion collisions provide a strong source of equivalent pho- 
tons up to very high energies. This offers the unique possibility to study 
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photon-hadron (nucleus) and photon-photon processes in hitherto inaccessi- 
ble regions, see [2]. Many UPC processes like GDR excitation are now recog- 
nized as useful for practical matters like luminosity measurement and impact- 
parameter dependent triggers. In this case, as was already shown explicitly in 
[46] one emphasizes the hard part of the spectrum, due to the selection of low 
b values. 
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